Analyzing the ecological relations of technology innovation of the Chinese high-tech industry based on the Lotka-Volterra model

Technology innovation has become an important driving force of economic and social development and has received wide attention from academics. Most scholars mainly take technology innovation as an overall variable to explore its impact on the economy and society. The main contribution of this study is to open the black box of technology innovation and introduce the lotka-Volterra model to explore the internal structure of technology innovation in the Chinese high-tech industry and to analyze the ecological relationships, evolutionary trends, equilibrium states of six technology innovation species including independent innovation (II), technology import (TI), research & development (RD), technology renovation (TR), foreign technology acquisition (FTA), and domestic technology purchase (DTP). The results of the study show that, First, the ecological relationship between prey and predator is observed between RD and TR, DTP and FTA, and II and TI. Second, no equilibrium state is observed between TD and TF and II and TI. Third, an unstable equilibrium state is observed between RD and TR.


Introduction
In the context of the COVID-19 pandemic, technology innovation plays a very important role in economic and social development as well as coping with global emergencies [1]. How to carry out technology innovation more effectively has attracted extensive attention from scholars [2]. High-tech industries, which are knowledge and technology-intensive, have long been the main carrier and driving force of technology innovation in China [3]. Based on the subjective desire for industrial upgrading or the objective need of the market, technology innovation has become an important driving force for the transformation and upgrading of high-tech industries and the enhancement of their core competitiveness [4]. The sources of technology innovation in high-tech industries include both internal and external aspects [5]. Internal technology sources include independent research and development (RD) and technological renovation (TR) [6,7] (FTA) and domestic technology purchase (DTP) [8,9]. The ecological relationships of technology innovation among RD, TR, FT, and DT not only affect the efficiency of technology innovation output of high-tech industries, but also have a great impact on the innovation ability, transformation, and upgrading of high-tech industries. Correctly analyzing and comparing the ecological relationships of the technology innovation in China's high-tech industry has become a real challenge that hinders the healthy and sustainable development of this industry.
Most of the existing studies have focused on analyzing technology innovation as a black box-like variable from two perspectives. On the one hand, some scholars have focused more on the role of technology innovation in the economic system. Jian et al. (2010) established the interaction between network relationships, trust, knowledge sharing, and technology innovation performance based on structural equation modeling [10]. Jiang et al. (2012) argued that high-quality technology innovation can promote technology standards and the widespread implementation of high-level technology standards can promote technology innovation [11]. Iranmanesh (2014) discussed the antecedents and outcomes of green technology innovation adoption in Malaysian transport firms [12]. Verbano and Crema (2016) investigated the role of technology innovation strategies on intellectual capital development and hence on innovation performance of manufacturing SMEs [13]. Feng and Yuan (2016) analyzed the impact of technology innovation and spillovers on the carbon intensity of human well-being using panel data of 30 provinces in China from 2005 to 2010 [14]. On the other hand, other scholars have concentrated on the factors for technology innovation. Shi and Zhu (2014) empirically examined the relationship between age, political connections, and technology innovation outputs of the pharmaceutical industry in China [15]. Sun (2015) explored the role of government in the technology innovation process in China from the perspective of strategic entrepreneurship [16]. Guo et al. (2016) argued that four types of factors, classified as innovation policy, innovation input, innovation capacity, and innovation organization, may influence the transformation of coal resource-based economies [17].
However, as a complex concept, technology innovation consists of several components depending on the source of the technology. In recent years, studies have discussed the components of technology innovation and their interrelationships. Kim and Stewart (1993) examined the relationship between domestic R&D and technology imports [18]. Lee (1996) analyzed the relationship between technology imports and the R&D efforts of Korean manufacturing firms through a two-stage approach [19]. Yu et al. (2016) conducted a study on the relationship between technical sources and innovation output of the Chinese high-tech industry [20]. But, the study of the relationship between different types of technology innovation did not cover all types, and the internal organic structure of technology innovation deserves further exploration.
Therefore, from an analytical point of view, opening the black box of technology innovation and exploring the internal organic structure of technology innovation is conducive to analyzing the internal connections and interactions between various types of technology innovation. In terms of research methodology, the static, one-way evaluation method of the relationship between factors related to technology innovation cannot conform to the dynamic development trend presented by technology innovation [21]. It is necessary to conduct a dynamic, ecological analysis of the two-way relationships of technology innovation. The Lotka-Volterra model is the most commonly used model for analyzing population competition and is very suitable for analyzing the evolution of ecological relationships of technology innovation populations.
Based on the above analysis, this study further opens the black box of technology innovation and divides it into two parts: independent innovation (II), including research & development (RD) and technology renovation (TR); and technology imports (TIs), including foreign technology acquisition (FTA) and domestic technology purchase (DTP), explores the complex dynamic relationship between the four types of technology innovation, and constructs a theoretical model among the four populations competition. We introduce the Lotka-Volterra model to explore the complex ecological relationships among technology innovation types in Chinese high-tech industries based on the data of technology innovation in China from 1995 to 2015. This study provides a theoretical basis for optimizing the allocation of technology innovation resources in high-tech industries and improving the efficiency of technology innovation. The key questions to be addressed in this study are as follows. What is the internal structure of technology innovation? What is the ecological relationship between the four types of technology innovation? What strategies should governments and enterprises adopt to strengthen technology innovation?
The remainder of this paper is organized as follows. Section 2 provides a short review of ecological relationships and the Lotka-Volterra model, and introduces the empirical model based on the Lotka-Volterra model and describes the dataset; Section 3 provides the empirical estimates and discusses the results; Section 4 presents the concluding remarks, policy implications, and limitations of this study briefly.

Ecological relationships and Lotka-Volterra model
Based on the theory of innovation ecosystems [22,23], the six types of technology innovations are regarded as six species that constitute a technology innovation ecosystem. We adopt the definition of ecological relationships (ERs) to discuss the cooperative and competitive interactions among the technology innovation species. Therefore, it is crucial to understand and consider the ecological relationships [24] that influence the balance and dynamics of technology innovation ecosystems [25].
Currently, descriptive mathematical models have been developed to predict the spread of species innovations and population development, such as the Gompertz model [9], logistic model [26], Bass model [27], GM (1, 1), and Lotka-Volterra model [28]; however, only the Lotka-Volterra model can examine the ERs of species.
The Lotka-Volterra model, first proposed by the American ecologist Lotka and the Italian mathematician Volterra [29], was derived to explore the growth curve-based interactions between two competing species, and it is often applied to describe the growth and decline of two competing groups in an ecosystem [30]. In recent years, this model has been applied in the field of socioeconomics to examine the issues of economic growth and population control [31]. Modis (1997) applied the model to study the dynamic relationship between competitors in a finite space [32]. Fu et al. (2017) used the Lotka-Volterra model to study the relationship between rural and urban Internet users in China [33]. Chakrabarti (2016) applied the Lotka-Volterra equations to model the evolution of the technology frontier [34]. Cordes and Schwesinger (2014) used an extended Lotka-Volterra model to capture the competitive interactions between two technologies [35]. Zhang and Lam (2013) introduced the Lotka-Volterra model to study the interaction between maritime sectors [36]. Lin (2013) applied the Lotka-Volterra model to analyze the competitive relationship between mobile cellular broadband and fixed broadband [37]. This study assumes that the technology innovation ecosystem conforms to the fundamental Lotka-Volterra model and that the ecological relationships can be described by corresponding mathematical equations.

Model analysis
According to the literature [31,37,38], the ecological relationships between two varieties of technology innovation in China's high-tech industry can be described by the Lotka-Volterra model using Eqs (1) and (2): In Eqs (1) and (2), X 1 and X 2 represent two technology innovation species; dX 1 (t)/dt and dX 2 (t)/dt denote the annual scale of X 1 and X 2 at each year t, respectively; X 1 (t) and X 2 (t) are the cumulative scales of X 1 and X 2 up to year t, respectively; X 1 (t)X 1 (t) and X 2 (t)X 2 (t) refer to the same technology innovation species interacting with itself; X 1 (t)X 2 (t) and X 2 (t)X 1 (t) denote two technology innovation species interacting with each other; α ij is used to determine the competition and cooperation of two technology innovation species; α 10 and α 20 are the logistic parameters representing the intrinsic growth rates of X 1 (t) and X 2 (t), respectively; α 11 and α 22 are the limitation parameters of the niche capacity for X 1 (t) and X 2 (t), respectively; and considering the possible endogeneity between variables, α 12 and α 21 are the interaction parameters representing how two technology innovation species affect each other. Modis (1999) summarized the possible types of ecological relationships of two species according to the signs of parameters α 12 and α 21 [39], as shown in Table 1.

Discretization
Since Eqs (1) and (2) are continuous-time models and the technology innovation scale in this study uses discrete data, the continuous Lotka-Volterra model must be converted into a discrete-time version. Currently, there are three discretization methods for the Lotka-Volterra model: the Leslie method, the log-integral method, and the gray method. In this study, the Leslie method was used for discretization, and the other two methods (See S1 Appendix for the specific formula derivation) were subsequently used for robustness testing.
According to the Leslie method [40], Eqs (1) and (2) can be transformed into discrete Eqs (3) and (4): In Eqs (3) and (4), b ii and b i0 are logistic parameters when only one species i lives alone and b 12 and b 21 refer to the degree of the effect of the relationship between two species on their respective growth rates. The relationship between the parameters in Eqs (1) and (2) and those in Eqs (3) and (4) are denoted as Eq (5). Note that the sign of α ij is the same as that of b ij : To perform regression calculations, Eqs (3) and (4)) can be transformed into Eqs (6) and (7): Let Eqs (6) and (7): can be converted into Eqs (9) and (10): Using the linear least square method, we can estimate the parameters in Eqs (8) and (10).

Model prediction ability calculation
This study uses the percentage error (PE), mean absolute percentage error (MAPE), mean absolute error (MAE), and root mean square error (RMSE) to measure the prediction ability of the model. PE and MAPE, MAE, and RMSE can be calculated using Eqs (11)~ (14): RMSE ¼ ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi 1 n In Eqs (11)~ (14), n represents the number of observations, A t is the actual value, and P t is the predicted value of A t . For all the above indexes, a smaller value implies a higher prediction accuracy. According to Lewis [41], the forecast capability levels of the MAPE are shown in Table 2.

Equilibrium analysis
Ecological relationships of technology innovation will reach a stable equilibrium when the volume of each species becomes constant, which means that the results of Eqs (1) and (2) are both zero [42][43][44][45]. We can obtain the equilibrium Eqs (15) and (16): (15) and (16) can be transformed into Eqs (17) and (18): In Eqs (17) and (18), X 1 is the horizontal coordinate and X 2 is the vertical coordinate. The evolutionary trend graph of X 1 and X 2 can be drawn in the plane rectangular coordinate system, according to which the evolutionary trends of the ecological relationships of technology innovations can be analyzed.
By solving Eqs (17) and (18), we can obtain the equilibrium points, as shown in Eq (19): a 11 a 20 À a 10 a 21 a 12 a 21 À a 11 a 22 However, not every equilibrium point is stable under the condition of environmental change. We applied the Jacobian matrix to calculate the Jacobian eigenvalue of each group and assess which equilibrium point of each group is stable. The condition for the stability of an equilibrium point is that both real numbers of eigenvalues are negative [34]. The Jacobian matrix is expressed in Eq (20):

Variables and data
Since the Lotka-Volterra model is more suitable for representing the relationship between two variables and two types of technology innovation are observed in the high-tech industries, namely, II, which can further be divided into RD and TR, and TI, which can further be divided into DTP and FTA, this study investigates the ecological relationships of three pairs of technology innovation: RD and TR (denoted ER1), DTP and FTA (denoted ER2), and II and TI (denoted ER3). The research framework of the ER of technology innovation is shown in Fig 1. The names, codes, units, definitions, and descriptions of all variables are shown in Table 3.
The original data in this paper are obtained from the China Statistical Yearbook on the hightech industry (1995~2015). Since technology innovation is a long process and requires a long time to accumulate, considering the availability of data, this study takes 2015 as the cut-off point and uses technology innovation data within 20 years as the sample. Considering the effect of the annual price level of technology innovation expenditure, the original data were deflated using the annual CPI index in 1995, and the data were obtained from the China Statistical Yearbook (1995~2015). Due to the data requirements of the Lotka-Volterra model, we also calculate the cumulative data for six variables of technology innovation in high-tech industries. The data of technology innovation in high-tech industries in China from 1995 to 2015 are shown in Table 4.

Results and discussion
In this section, we first analyze the models proposed in Sec. 2 with the accumulated expenditure data for technology innovation of China's high-tech industry (1995~2015), Then, we use

PLOS ONE
The ecological relations of technology innovation of the Chinese high-tech industry the criteria proposed in Sec. 3.7, i.e., MAPE, MAE, and RMSE, to compare the prediction accuracy of the three discretization methods.

Correlation analysis
The correlation analysis results for the six technology innovations in Table 5 show that the correlation coefficients of the six variables are very close to 1 and significant at the 0.01 level, implying that there is a clear correlation between the six variables, thus providing a foundation for studying the ecological relationships between variables.

Regression analysis
According to the data in Table 4, we estimate the parameters of the Lotka-Volterra model using three discretization methods, namely, the Leslie method, the log-integral method, and the gray method. In the analysis process of ER1, RD was designated X 1 and RD was designated X 2 . In the same way, DTP was designated X 1 and FTA was designated X 2 for ER2 and II was designated X 1 and TI was designated X 2 for ER3.
In the regression analysis of the Lotka-Volterra model with the Leslie method, we first estimate the parameter C ij in Eqs (9) and (10) with the data in Table 4, substitute the values of C ij into Eq (8) to obtain parameter b ij of Eqs (3) and (4), and substitute the value of b ij into Eq (5) to calculate parameter α ij of Eqs (1) and (2). Table 6 shows the estimation and calculation results of each parameter of the Lotka-Volterra model with the Leslie method, and the findings show that the statistical significance of the estimated parameters is generally remained below 1%, with only a few parameters remaining below 5% according to the t-statistics of parameter C ij . The model fit is approximately 0.7 according to the R 2 value, and the limit of F is below 0.01 according to the F value.
The parameters and model test results in Table 6 show that the parameters are significant according to the p-value of the t-test, the model has an extremely good fitting degree according to the values of R 2 and the regression equation has obvious linear significance according to the p-value of the F-test.

Ecological relationship analysis
Eqs (1) and (2) are based on the Lotka-Volterra model and contain all of the fundamental parameters that affect the growth rates of the six technology innovation scales. The ecological relationships among the six technology innovations can be determined through the parameter α ij based on the intrinsic growth rate, niche capacity limitations, and interaction, as shown in Fig 2. (1) Intrinsic growth rate. According to parameters α 10 and α 20 in Table 6, the intrinsic growth rate of six technology innovation species can be analyzed. The intrinsic growth rate parameters α 10 and α 20 of the six technology innovation species are all positive, indicating that the scale of these innovations is consistently growing from 1995 to 2015. In addition, the parameter α 10 for RD is 0.6347199313, which is the maximum value among the six intrinsic growth rate parameters, thus implying that the average growth rate of RD is the fastest, while

PLOS ONE
The ecological relations of technology innovation of the Chinese high-tech industry the parameter α 10 of DTP is 0.2910435054, which is the minimum value among the six intrinsic growth rate parameters, thus implying that the average growth rate of DTP is the slowest. The results of the analysis are consistent with the actual data in Table 4, validating the applicability of the Lotka-Volterra model for analyzing technology innovation.
(2) Niche capacity limitations. According to parameters α 11 and α 22 in Table 6, the niche capacity limitations of the six technology innovation species can be examined. The values of Parameter α 11 for RD, DTP and II are equal to 0.0000983798, 0.0018554485, and 0.0000244737, respectively, which are all positive. The regression results indicate that the niche capacity has a promoting effect on RD, DTP, and II. In addition, the value of DTP is the maximum and that for II is the smallest. Parameter α 22 values for TR, FTA and TI are equal to -0.0003054643, -0.0004383309, and -0.0004122655, respectively, which are all negative. The regression results demonstrate that the niche capacity has an inhibitory effect on TR, FTA, and TI. The value of FTA is the maximum and that for TR is the smallest. These analytical results may be interpreted as follows. Although RD, DTP, and II, are the core power and significant for the sustainable development of the Chinese high-tech industry, the scale of these innovations is not sufficient to forcefully transform and quickly upgrade the Chinese high-tech industry. For the technology innovation system, niche capacity has a promoting effect on RD, DTP, and II. However, TR, FTA, and TI are not core power of the Chinese high-tech industry and represent supplements for RD, DTP, and II, respectively. The scales of TR, FTA, and TI are large enough to supplement RD, DTP, and II. Under the context of insufficient RD, DTP, and II and the systematic view of technology innovation, it is reasonable that niche capacity will have an inhibitory effect on TR, FTA, and TI.
(3) Interaction. Based on the regression results, the interaction parameters α 12 and α 21 in Table 6 and the type of the ecological relationships in Table 1, it is possible to examine the interaction between the six technology innovation species. For RD and TR (ER1), the parameter α 12 is -0.0004267798, which is negative for RD, meaning that TR inhibits the subsequent scale of RD, while parameter α 21 is 0.0000711148, which is positive for TR, implying that RD promotes the subsequent scale of TR. These analytical results demonstrate that RD is a prey species and TR is a predator species, with RD serving as direct food for TR.
For DTP and FTA (ER2), the parameter α 12 is -0.0003425330, which is negative for DTP, implying that FTA inhibits the subsequent scale of DTP, while the parameter α 21 is 0.0011854446 and positive for DTP, meaning that FTA promotes the subsequent scale of DTP. These analytical results demonstrate that DTP is a prey species and FTA is a predator species, with DTP serving as direct food for FTA.
For II and TI (ER3), the parameter α 12 is -0.0003428341, which is negative for II, implying that TI inhibits the subsequent scale of II, while the parameter α 21 is 0.0000248746, which is positive for TI, implying that II promotes the subsequent scale of TI. These analytical results demonstrate that II is a prey species and TI is a predator species, with II serving as a direct food for TI.
The results of the above analysis indicate three conclusions. First, all six technology innovations increased from 1995 to 2015. Second, RD, DTP, and II promote themselves while TR, FTA, and TI inhibit themselves. Third, the ecological relationships between RD and TR, DTP and FTA, and II and TI are all prey-predator relationships. Therefore, it can be seen that the population size change of each species of technology innovation is not only influenced by itself, but also by the number of other species, and the study of the relationship between the components of technology innovation from the ecological perspective is consistent with the need for technology innovation.

Evolutionary trend analysis
For RD and TR (ER1), by substituting parameter α ij of RD in Table 6 into Eq (17), we can obtain Eq (21), and by substituting the parameter α ij of TR in Table 6 into Eq (18), we can obtain Eq (22), L 1 : 0:6347199313 þ 0:0000983798RDðtÞ À 0:0004267798TRðtÞ ¼ 0; ð21Þ L 2 : 0:3941010829 þ 0:0000711148RDðtÞ À 0:0003054643TRðtÞ ¼ 0: Eqs (21) and (22) Table 4. Since both RD and TR are greater than zero, only the region of the first quadrant is of practical interest for the analysis of the evolutionary trend of RD and TR. L 1 and L 2 intersect in the first quadrant, indicating an equilibrium point for RD and TR. In addition, L 1 and L 2 divide the first quadrant into four regions. For any point in region I, we can obtain the following inequalities (23) In Eqs (23) and (24) show that, dRD(t)/dt<0 and dTR(t)/dt<0 according to Eqs (15) and (16). Therefore, in region I, RD and TR will decrease. Similarly, in region II, dRD(t)/dt>0 and dTR(t)/dt<0, indicating that RD will increase while TR will decrease. In region III, dRD(t)/ dt>0 and dTR(t)/dt<0, indicating that RD and TR will increase. To the right of the equilibrium point, to the left of L1, and the right of L2 is region IV, where dRD(t)/dt<0 and dTR(t)/ dt>0, indicating that RD will decrease and TR will increase.
The results of the above analysis indicate three conclusions. First, an equilibrium point emerges between RD and TR, but its stability must be verified, while no equilibrium point emerges between DTP and FTA, II and TI. Second, RD, DTP, and II will decrease, while TR, FTA, and TI will increase. Thus, it can be seen that opening the black box of technology innovation and considering each type of technology innovation as a species can provide a clearer study of the evolutionary trends of the components of technology innovation and provide worthy ideas for technology innovation research.

Equilibrium stability analysis
By substituting the values of the parameter α ij in models ER1, ER2, and ER3 into Eq (19), we can obtain the equilibrium points of ER, as shown in Table 7. The actual economic data show that the data must be positive and X 1 >0, X 2 >0, according to which we can confirm that only

PLOS ONE
ER1 has a realistic equilibrium point between RD and TR while ER2 and ER3 have no realistic equilibrium point, thus verifying the results in section 6.
Next, the stability of the equilibrium point of the ER1 model is determined based on the eigenvalues of the Jacobian matrix. Substituting the values of the parameter α ij and the equilibrium point of model ER1 into Eq (19), we can obtain the eigenvalues of the Jacobian matrix of model ER1, i.e., 1.77 and 0.19, according to which the real part of these eigenvalues is greater than zero. The eigenvalue does not satisfy the stability condition, thus implying that the stability of the equilibrium point of ER1 is unstable. Based on the above equilibrium analysis, we find that no equilibrium point appears between TD and TF and II and TI, and an unstable equilibrium point appears between RD and TR.

Robustness test
In this study, the regression analysis of the model was done again using the log-integral method and the gray method to verify the robustness of the results. Table 7 shows the results of the estimation and calculation of each parameter of the Lotka-Volterra model with the logintegral method, and the results of the study imply that the statistical significance of the estimated parameters is generally staying below 1%, a few parameters stay below 5%, and only one parameter stays above 10% according to the t-statistic of the parameter C ij . Based on the R 2 value, the model fit is approximately 0.8, and the limit of F value is below 0.001. Table 8 shows the results of the estimation and calculation of each parameter of the Lotka-Volterra model with the gray method, and the results show that the statistical significance of the estimated parameters mostly stays below 1% level according to the t-statistic of the parameters C ij, a few parameters stay below 5% or 10%, and only one parameter stays above 10%. The model fit is approximately 0.8 according to the R 2 value, and the bounds of F are below 0.001 according to the F value.
The results of the parameters and model tests in Tables 8 and 9 show that the parameters are significant according to the p-value of the t-test, that the model fits very well according to the value of R 2 , and that the regression equation is significantly linear according to the p-value of the F-test. Thus, the results show a high degree of robustness.

Comparison of the forecast ability levels
Substituting the actual data in Table 4 and parameter values of b ij in Tables (6) to (8) obtained by the three discretization methods into Eqs (3) and (4), we can calculate the estimated values of the six technological innovations by the three discretization methods. Then, substituting the estimated values and the actual data into Eqs (11)~(14), we can obtain the values of MAPE, MAE, and RMSE for the six technological innovations by the three discretization methods, as shown in Table 10. Among the three discretization methods, the Leslie method has the smallest MAPE value and the same MAE and RMSE values for the six technological innovations, which means that the Leslie method is the most suitable for the Lotka-Volterra model among the three discretization methods. Therefore, this paper analyzes the ecological relationships, evolutionary trends, and equilibrium states of the six technological innovations using the estimated parameters of the Lotka-Volterra model with the Leslie method.

Conclusions
This study investigates the ecological relationships among technology innovations in China's high-tech industry. Based on the Lotka-Volterra model and the data from 1995-2015 in the Statistical Yearbook of China's High-Tech Industries, we derive interesting and important conclusions. First, the Leslie method is the best discretization method for the Lotka-Volterra model and the Lotka-Volterra model is suitable for describing and predicting the scale changes of six technology innovation species in Chinese high-tech industries. Second, the ecological relationships between RD and TR, DTP and FTA, II and TI are prey-predator relationships, and RD and DTP and II will be greatly reduced, while TR and TF and TI will be gradually increased. Third, there is no equilibrium point between DTP and FTA, II and TI, and an unstable equilibrium point between RD and TR.

Implications
Based on the above analysis, this study draws the following three policy implications: 1. The government should strengthen its support and subsidies for high-tech industries. While reducing R&D uncertainty and risk, the government should also focus on the implementation of anti-monopoly rules and strategies to prevent large downstream buyers from abusing market dominance and protect the innovation resources of enterprises, and create a good market competition environment for enterprise technology innovation.
2. The government should continue to promote globalization and introduce R&D and innovation foreign direct investment under local conditions. While improving the intellectual property protection system, the government should introduce high-quality foreign investment in a targeted manner, give full play to the role of foreign technology innovation in promoting the technological progress of domestic enterprises, and reasonably optimize the level of foreign-funded competition in different industries rationally to provide nourishment for improving local independent innovation capability.
3. High-tech enterprises should implement technology diversification strategy, and realize both collaborative and exploratory innovation. Enterprises should make full use of the advantageous knowledge resources of the members of the collaborative innovation network, give full play to the important organizational capabilities of product innovation, optimize the structure through technological diversification, and open up various links of resource integration.

Limitations
Two limitations should be taken into consideration in subsequent research. First of all, this study only examined the input data of technological innovation because it is difficult to obtain the output data of technological innovation, which can also indicate the development trend of various technological innovations. Moreover, we ignore the life periodicity of technology innovation of the high-tech industry. It would be meaningful to study the evolution of ecological relationships among different varieties of technology innovation during various periods of the life cycle.
Supporting information S1 Data.